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A general theory is derived for the moments of the first passage time of a one-dimensional Markov 
process in presence of a weak time-dependent forcing. The linear corrections to the moments can 
be expressed by quadratures of the potential and of the time-dependent probability density of 
the unperturbed system or equivalently by its Laplace transform. If none of the latter functions 
is known, the derived formulas may still be useful for specific cases including a slow driving or 
a driving with power at only small or large times. In the second part of the paper, explicite 
expressions for mean and variance of the first passage time are derived for the cases of a linear or 
a parabolic potential and an exponentially decaying driving force. The analytical results are found 
to be in excellent agreement with computer simulations of the respective first-passage processes. 
The particular examples furthermore demonstrate that already the effect of a simple exponential 
driving can be fairly involved implying a nontrivial nonmonotonous behavior of mean and variance 
. as functions of the driving's time scale. 

Keywords: first-passage time; driven stochastic systems 



d . I. INTRODUCTION 

C/3 , 

■ One of the key results in the theory of stochastic processes are the quadrature expressions for the moments 
of the first passage time (FPT) in case of a one-dimensional Markov process 0, • One of the assumptions 
made in the seminal papers by Pontryagin et al. [l| and Siegert 0] is the temporal homogeneity of the process: 
except for the driving Gaussian white noise the system is subject to a force without time dependence. In many 
instances this assumption does not hold and an extension of the classic theory to the case of a time-dependent 

Q ' force is desirable. 

O , In the context of many noise-induced phenomena, for instance, the presence of an additional deterministic or 
stochastic perturbation is essential and its effect on the passage time statistics is of foremost interest. In case 
y—i \ of resonant activation (RA) jj, 0, la, E|, the additional driving is a stochastic process with state-dependent 

J> • amplitude. For the problems of stochastic resonance (SR) (see Ref. and references therein) and coherent 
stochastic resonance (csr) HUGS lllj the driving is commonly a deterministic periodic signal. The key 
feature of both RA and CSR is a minimum in the mean FPT or mean exit time attained at a finite "optimal" 
value of the forcing's time scale (e.g. correlation time or driving period). Similarly, SR is realized if a time-scale 
matching condition between the forcing period and the FPT of the unperturbed system (an inverse Kramers 

■ rate) is met. 

My own interest in the general problem originates in studies of stochastic neuron models involving an expo- 
nentially decaying time-dependent forcing. For these models, the FPT corresponds to the interspike intervals 
(ISIs) separating the neural discharges (action potentials or spikes) by means of which neurons transmit and 
process signals. The exponentially decaying perturbation in these models arises from slow ionic currents that 
are driven by the spiking activity of the neuron itself [l^ . An exponential perturbation is also obtained via a 
simple transformation of models with a decaying threshold; such models have been employed to reproduce the 
statistics of certain sensory neurons 0, ^| . In general, the effect of an exponentially decaying driving on the 
£j \ FPT statistics is poorly understood up to now in contrast to the frequently studied case of periodically forced 
stochastic neuron models (see e.g. Refs. [H[ll [H[ll [H E3). 

In order to study the change of the FPT statistics induced by a time-dependent driving, researchers have 
used two different analytical approaches. First, keeping the potential shape, the driving force as well as 
the boundary and initial conditions as simple as possible allows in some cases for an exact, in others for an 
approximate solution. This kind of approach was pursued in the late 1980's H, 0, ^ and during the 1990's 
0,13, El ITU ITtH (see also Ref. j2^| , ch. 4 for further examples) ; in most cases (piecewise) linear potentials and a 
periodic or dichotomous driving were considered. Later several researchers proposed semi-analytic procedures 
to solve these FPT problems EE 26] . It should also be noted that there exists a considerable mathematical 
literature on the FPT problem with time-dependent deterministic forcing, most of which is devoted to specific 
simple cases that are exactly solvable by a transformation to a time-homogeneous system (see Refs. p7ll28ll29| 
and further references cited in these papers). 

The other analytical approach consists of a weak noise analysis of driven first-passa ge p rocesses. Here the 
only assumption made about the potential is the existence of a metastable state. Refs. |3(J, IU |32|, y3( focused 
on the escape rate out of a metastable one-dimensional potential (the inverse mean first passage time for a 
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quasi-equilibrium initial condition) in the small-noise limit and aimed thus at a generalization of the famous 
Kramers problem 0, |3j| for periodic [3(], On], or stochastic 0, @ driving. Although both approaches led 
to significant progress, there are also many instances where one wishes to go beyond the limit set by assuming 
weak noise or piecewise linear potentials and where also higher moments and not only the mean FPT (or its 
inverse) may be of interest. 

Here I give an extension of the classic concept due to Pontryagin pj and Siegert Q of calculating the FPT 
moments for a general potential to the case that includes a weak time-dependent driving. I focus on a 
deterministic driving function; the approach might be, however, also helpful in situations with an additional 
stochastic driving. I will use a perturbation calculation that requires a weak driving and will consider the first 
(linear) correction term for each moment. No restrictions apply with respect to the noise intensity or the time 
scale of the driving, although the range of validity of the theory, i.e. the question how "weak" the driving has 
to be, may also depend on these parameters. My most general result relates the linear corrections to each 
moment of the FPT to quadratures of the time-dependent probability density Po(x,t) of the unperturbed 
system or to quadratures of the Laplace transform of this function. An alternative formulation provides the 
corrections in terms of an infinite sum of quadratures of known functions. These results cannot be applied in 
the general case (arbitrary driving function and arbitrary potential shape), since Po(x,t) and/or its Fourier 
transform are not known for most potentials and the numerical evaluation of an infinite sum of quadratures 
is in general not feasible either. For many important cases, however, including a slow driving, a driving with 
power at only small times or a driving with power at only large times, the derived general formulas can be 
useful for the FPT problem in a potential of arbitrary shape. Furthermore, I will show that for two specific 
potential shapes (linear and parabolic) and an exponentially decaying driving force, the derived theory yields 
valid (and in part strikingly simple) results for the FPT's mean and variance for arbitrary time scale of the 
driving (slow, moderate, or fast compared to the mean FPT of the unperturbed system). The explicite results 
derived for an exponential driving function can be applied to the neurobiological problems mentioned above; 
this will be done elsewhere. Since the exponential decay of drift parameters is frequently encountered in many 
situations, the results may be also helpful for the study of other physical systems. 

This paper is organized as follows. Sec. [H] starts with an introduction of the problem and of the quantities 
of interest. Next I present an derivation of alternative quadrature formulas for the FPT moments in the 
time-homogeneous case. These quadratures are entirely equivalent to the classic formulas by Siegert Q as 
shown in the appendix. The alternative approach will be the basis for the perturbation calculation of the 
time-inhomogeneous first-passage problem leading to the general relation between the corrections of the 
moments and the quadratures of the unperturbed probability density. In sec. IIIII I derive explicite results 
for a linear and a parabolic potential and an exponential driving force. These analytical results will be 
compared to simulations of the two systems in sec. IIVI Here I will show that, remarkably, similar to the case 
of periodic driving, a nontrivial behavior of the FPT's mean and variance with respect to the time scale of the 
driving (i.e. the decay rate of the exponential driving) is possible. The mechanisms for these "resonances" 
will be discussed. Sec^summarizes the findings and discusses further applications and extensions of the theory. 



II. GENERAL THEORY 

A. Langevin and Fokker-Planck equations 

Starting point of my consideration is the Langevin equation for a one-dimensional escape process given by a 
potential U(x), a white-noise driving of intensity D, and an additional weak time-dependent forcing es(t): 

x = -U'{x) + es{t) + V2D£{t) (1) 

Without loss of generality I consider an initial value at zero (x(0) = 0) and ask for the first-passage time to a 
point xe to the right of the origin (xe > 0). The only restriction for the potential is that I exclude potentials 
that allow for an escape toward minus infinity. Furthermore, I do not specify the forcing s(t) but take for 
granted that its effect on the FPT statistics is weak by virtue of the small parameter e. 

Instead of using the Kolmogorov (backward) equation as it is commonly done in the treatment of FPT problems 
0, , I shall employ the Fokker-Planck (forward) equation (FPE) governing the probability density of x 



d t P(x,t) = d x [(U'(x) - Es(t))P(x,t) + Dd x P(x,t)]. 



(2) 
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The FPT problem stated above determines the initial condition (at t = the variable x is with certainty at 
x = 0) and the boundary condition (absorbing boundary condition at x = xe) 

P(x,0) = S(x), (3) 
P(x E ,t) = 0. (4) 

It is well known that there is a simple relation between the FPT density and a quantity derived from the 
probability density P(x, t): the FPT density is given by the time-dependent probability current at the absorbing 
boundary (see, for instance, Ref. |37jl 

$(T) = J(x E> T) = -[(U'(x)-es(T))P(x,T) + Dd x P(x,T)}\ x = XB 

= —Dd x P(x, T)\ X=XE . (5) 

Provided Q(t) is known, one can calculate the moments of the FPT by 

oo oo 

(T n (0^x E )) = J dT T"$(T) = -D J dT T n d x P(x,T)\ x=XE . (6) 
o o 

For certain problems it may be desirable to know also the Laplace transform of $(T). This function can be 
expressed by the Laplace transform of the probability density p(x,p) as follows 

oo 

pip) = Jdt e-f*$(t) = -Dd x p(x,p) . (7) 

X = XE 

By means of the Laplace transform the moments can be calculated as follows 

(8) 



(T n (0^x E )) = (-ir^ P (p) 



p=0 



However, even in the absence of a time-dependent driving (e — 0), solving directly for one of the functions P(x, t), 
$(T), or p(p) is possible only in a few simple cases including constant and linear potential shapes. Nevertheless, 
for calculating the moments of the FPT we are not restricted to use eq. (0 or eq. JHJ- Remarkably if e = 0, the 
moments of the FPT can be calculated exactly for an arbitrary potential shape U(x). The general formula for 
the n-th moment of the passage from a general initial point x = a to an absorbing boundary at x = b is given 

by i 



OT'Ci - - /,))„ = g [ dye u M' D f dxe~ u ^ ' D {T^ 1 {x -► 6)> (9) 



oo 



where the index "0" indicates that e — (later on, this convention will be also applied to the functions 
P(x,t),p(x,p) and pip))- In order to calculate the n-th moment one has to solve for the lower moments as 
functions of the initial point first. The hierarchy of quadratures is completed by stating that for obvious reasons 
(T°(o-.6)) = (l) = l ; 

The aim of this paper is to extend these expressions to the case of a weak time-dependent driving function es(t). 
In other words, I seek for the linear correction term to the n-th moment denoted j n {xE) such that 

(T") = (T") + ej n (xE)- (10) 

In particular, once the corrections to the first two moments have been calculated, mean and variance of the 
FPT will be to linear order in e given by 

(T) - Ma + ehixE), (11) 
{AT 2 ) = {T 2 ) - (T) 2 

= (AT 2 ) Q +e[j 2 (x E )-2(T)oj 1 (x E )}. (12) 

Later on, for specific systems, I will solely discuss these first two cumulants and the relative standard deviation 
(coefficient of variation) of the FPT, that is a function of them 



(13) 
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B. Moments of the first passage time in the autonomous case - the other way around 

Here I set e = 0. First I introduce the functions 



/>oo 

(3 n (x) = / dt t n P(x,t) (14) 
Jo 



and 



J n (x) = -(U'(x) + D-^)Pn(x) (15) 

On comparing eq. © and eqs. (|14f> and (|15|) it becomes apparent that j3 n and J n are related to the n-th moment 
of the FPT as follows 

(T") = -DpJx E ) = J n (x E ). (16) 

Here and in the following the prime denotes the derivative with respect to x. I now derive the general solutions 
for J n (x) that provides an alternative solution for the FPT moments by virtue of eq. (|16fl . 
Multiplying the FPE J2J with t n , integrating over t, and using integration by part on the l.h.s. of the equation, 
I obtain for the functions (i n {x) the following set of equations 

-S(x) = ±(U'(x)+D±)l3o(x), (17) 
-np n ^(x) = -^(U'(x)+D-^)0 n (x), (18) 
The boundary conditions can be inferred from those for P${x,t) 

0„(x E )=O and /^(-oo) = (with k = 0,1,2, • • •) (19) 
with P^\x) denotes the A:-th derivative. The solutions are straightforward 

X E 

(3 (x) = l e - u ^ D J dye u ^/ D e(y), (20) 

X 

xe y 

Pn{x) = ^ e - u ^/D J dye U(y)/D J dz ^ l{z)i n = l,2,--- (21) 

x -co 

where B(x) denotes the Heaviside jump function 38]. For the J n {x), I find 

-S(x) = -J' {x), (22) 
-n(3 n _ x {x) = -J' n {x), n = l,2,--- (23) 

The first equation leads immediately to a Heaviside function 

J (x) = Q(x), (24) 

To the second equation I apply the operator —(U'(x) + Dd/dx) which yields 

-nJ n - 1 {x) = U'{x)J l n {x)+Dj!l{x), n = l,2,-.- (25) 

From eqs. (|22|) . l|23(l . and l|19(l . I get boundary conditions for the J„ 

J' n {x E )=0, J„(-oo) = (26) 

The solution of eq. I)25|l obeying these conditions reads 

X xe 

Ux) = £ J dy e-^^ J dz e u M' D J n -i(z) (27) 
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This equation together with eq. (|24|l and eq. I|16|) provides an alternative, though not in the least simpler way 
of calculating the moments of the FPT. On comparing with eq. 10, I note the differences in the signs of the 
exponents, in the boundaries of integration, and in the first function of the hierarchy eq. Ij24(l (for eq. ©, the 
first of the functions was (T°(a — > b)) = 1). Nevertheless, the alternative quadrature expressions derived here 
are completely equivalent to eq. © as shown in the appendix. In particular for n — 1 eq. (|27|l and eq. (|24|l yield 



(T) = J\{x E ) = i 



° dx e - u ^ D I** dy e u ^/ D + 



30 JO 

XE dx e~ u ^/° r dy e u ^' D 



(28) 



Changing the order of integration yields then the same result as eq. (0 for n = 1, a = 0, b = x E - 

C. Including a weak time-dependent force 

I now turn to the case £ ^ 0. For the probability density P(x,t) obeying the FPE J2J, I make the following 
ansatz 

P{x,t) =P (x,t) + eq(x,t) (29) 

where the first term is the probability density for e = 0. Omitting all terms of order e 2 and higher, I find the 
following equation governing the second function 

d t q(x,t) = d x {U'(x)+Dd x )q(x,t) - s{t)d x P (x,t) (30) 

The boundary and initial conditions for this function can be inferred from those of P(x,t) and Po(x,t) 

P(x E ,t)=Po(xE,t) = -» q(x E ,t)=0 (31) 
P(x,0) =P o (x,0) = 5{x) -> q(x,0) =0 (32) 
P(-oo,t) = P (-oo,t) = q(-oo,t)=0 (33) 

Now I introduce the counterpart to the functions J n {x) corresponding to the perturbation q(x,t) 

poo 

j n (x) =-(U' + Dd x ) / dt t n q(x, t) (34) 



JO 

Knowledge of this function allows for the calculation of the n-th moment of the first passage time by the 
following formula 

(T n ) — (T") + ej n (x E )- (35) 

From eq. (SU, I find 



-j'Jx) = -n\ dt t n - 1 q{x,t)+ dt t n s(t)d x P (x,t) (36) 
Jo Jo 

From this equation and from eq. (|33(l and eq. (|31(l . I can conclude that 

i„(-oo) = (37) 

j' n (x E ) = dt t n s(t)d x Po(x E ,t) = -d x f n P (x,t) (38) 

Jo 

where I have used the integral operator T n that is defined by 

roc 

f n = dt s(t)t n (39) 
(function is multiplied by t n s(t) and then integrated). The solution for jo(x) is straightforward 

/>oo 

j {x) = -f P (x,t) = - dt s(t)P (x,t) (40) 
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The constant of integration in eq. i|40|) is zero because of the boundary conditions eq. (|31|l and eq. (|37|l . 
For n > 0, I apply the operator —(U' + Dd x ) to eq. 

U'j' n + Dj'l = -ni„_i - f„(L/' + Da x )P (x, t) = f n (x) 



(41) 



This equation has the same structure like those for the functions J n (x) of the unperturbed system. The difference 
lays in the inhomogeneities on the r.h.s. (abbreviated by f n (x)) and the different boundary condition for the 
derivative of j n at xe ■ 
The solution for the derivative is given by 



EM 

Jn(») = e D 



1 f EM 

"n + / dy e D f n (y) 



EM 
e d 



n I EM - EM 

I dy e d j„_i(y)-T n e d d y P {y,t) 



D 



U(x) 



dy e d dMy,t) 



U(x) 



dy e d d£P {y,t) 



(42) 



The terms in the last line cancel and the integration constant c n has to be chosen such that condition eq. (|38|l 
is met. I obtain 



j'Jx) = -f n d x P (x,t) + -e- u ^ D I dy e u W D 3n-li3l) 



(43) 



Another integration (for which the integration constant is determined by the boundary condition at — oo given 
in eq. (J23) yields 



j n (x) = -f n P (x,t) + ^ / dz e' u ^' D I dy e u ^l D j n - X {y) 



(44) 



This equation is the first important result of my paper. I recall that the corrections to the moments of the 
first passage time are obtained by taking j n at ig. Thus if the function P (x 1 t) is given, one may calculate 
the effect of the external driving on an arbitrary moment of the first passage time by a subsequent solution of 
all the lower moments. If P^{x,t) is not given (which is, unfortunately, usually the case) there are still several 
classes of tractable problems for which eq. (|4*4"}l is useful. These are discussed in the next subsection. 

D. Further simplification of the result for specific cases 

Eq. Q44f> involves integrals over the probability density of the unperturbed system multiplied with the time- 
dependent part of the drift and powers of t. Expanding the driving function in powers of t permits to express 
these integrals in terms of the known functions (3k from eq. I|21|l as follows 



f n P (x,t) = dt s(t)t n P {x,t) = ■ 

Jo fc=0 



,(fe) 



(0) 



Pk+n(x) 



by means of which I obtain 



jnix) 



D 



dz e- u ^/° 



dy e u ^ D Jn ^(y)-Y- 

k=0 



S W(0) 
k\ 



(45) 



(46) 



This formula is especially useful in case of a slow driving that can be for t ~ (T} described by just a few 
expansion terms . As can be expected, the zeroth term results in the static correction, what I briefly show 
now for the simplest case n = 1. Suppose a static driving s(t) = 1, then jo(x) = —f3o(x) and the linear correction 
to the mean FPT reads 



h, s taUc(x E ) = -i f dze- u ^ D fdye u ^ D My) 
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D 2 



dz e - u{z)/D / dx / dye u{v) Q(y) 



Z X 



xe y 

—ac 



(47) 



This correction is also obtained by considering the unperturbed system with a potential U(x) = U(x) — ex, 
writing down the mean FPT according to eq. JHJ with n = 1, a = 0, b = xe, and expanding the result up to first 
order in e. For a non-static but slow forcing, the first correction term describing a truly dynamical effect of the 
driving would be obtained by taking into account a finite s'(0), leading to a quadrature formula that involves 
Although the incorporation of higher non-static corrections is straightforward, note that the number of 
quadratures which have to be numerically solved increases by two with every term s^ n ^(0) that is taken into 
account. 

Besides a slow driving another important class of perturbations is given by drivings described by an exponential 
decay or a periodic function, both of which can be described by exp[— Xt] or a sum of such exponentials. In this 
case the term T n PQ(x,t) equals the n-th derivative of the Laplace transform po(x,p) of Po(x,t) with respect to 
the complex argument p taken at p = A 



T n P (x,t) 



dt e- pt t n P {x,t) 



jn 

= (-i) v^' p) 



(48) 



P =\ 



Using this form in eq. (|44|) is in particular of advantage if the function fto(x,p) is known but Po(x,t) is not 
known. For more general driving functions, this can be generalized as follows. Suppose the Laplace transform 
of the driving s(t) exists 



(w) = / dt e- pt s(t), (p>0) 



(49) 



Then it is possible to recast eqs. 11441 into the following form 



j n (x) = -f n po(x, -p) + - / dz e - u ^ D I dy e u ^l D j n -i{y) 



(50) 



where the operator T n is defined by 



1 f d n 
J~ n = - — / dp s{p)— — 
2ni J dp n 



(51) 



With a pure exponential driving, the term J r n po(x, —p) reduces to eq. (|48fl as can be shown by the calculus of 
residues. 

Further simplifications or approximations are possible by means of short-time or asymptotic approximations of 
Po(x,t) (see, for instance, Ref. |39|). if most of the driving power is at short or long times. Here I shall not 
further dwell on the general case but study the effect of a simple exponential driving on systems with linear or 
parabolic potential for which exact correction formulas for mean and variance can be found. 

III. THEORY FOR A SYSTEM WITH LINEAR OR PARABOLIC POTENTIAL AND 

EXPONENTIAL DRIVING 

In the following, I will focus on the corrections to the first two moments. These are given by 

3i{xe)=~ J dxe- u ^/ D ' dy e u ^ D f P (y,t) (52) 

— OO 
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h(x E ) = -~ I" dx e' u ^ D dy e u (yy° {f 1 P (y,t)+ 

i f daae"^^ f% 2 e ^ B f i%(» 2 ,t)} (53) 

U J — OO J X2 ) 

Furthermore, the following driving force is considered 

s(t) = Aexp[-At]. (54) 

Some remarks regarding this function are indicated. The function s(t) is normalized (integration over time 
yields one). There are two simple limits: (i) for A <C (To), the function tends to a static bias of amplitude A; (ii) 
for A — > oo the function approaches a S function that changes the initial position from to e. In both limits, 
the moments of the first passage time can be exactly calculated which gives us a mean to check the validity of 
the perturbation calculation. Furthermore, according to eq. the operators T n correspond to derivatives of 
the Laplace transform po(x, A) multiplied with A 

d n 

T n P (x,t) = A(-l)"— po(x,X) (55) 

Finally note, that the correction formulas can be looked upon as linear operations on the driving function. This 
implies that the correction to a driving consisting of a sum of exponentials equals the sum of the corrections to the 
single exponentials. In particular, this applies to a (possibly damped) cosine driving s(t) = exp[— (X + iuj)t] +c.c. 
with A, to € 3? and A > 0. 

In the following, I will furthermore use a parabolic potential 

x 2 

U(x) = b Y~ ax > ( b - °) ( 56 ) 

and will separately discuss the cases b — and 6^0. The former problem corresponds with a > to a biased 
random walk toward the absorbing boundary; there is no potential barrier present in this simple case and 
the first passage will take place in a limited time even at vanishing noise. In contrast to this, for 6^0 and 
a/b < xe there exists a metastable point (potential minimum of U(x)) to the left of the absorbing boundary; 
the first passage process is noise-assisted, i.e. for vanishing noise the passage time tends to infinity. I note that 
both cases are of particular importance in the neurobiological context, where the FPT corresponds to the so 
called interspike intervals generated by a perfect (b = 0) or leaky (b > 0) integrate-and-fire neuron stimulated 
by white noise |3^ |. 

Before I come to the specific cases, some further simplifications of the correction formulas eq. (|52|l and eq. (|53|l 
for the potential eq. (|56J) with arbitrary b will be performed. It will emerge that for a general parabolic 
potential the corrections can be entirely expressed by the Laplace transform po (A) of the FPT density of the 
unperturbed system instead by the function po^j A). 



A. Simplification of the correction formulas for arbitrary b 
Using eq. I|55|l with n = the correction to the mean FPT eq. I|52|l reads 

xe 

ji(x B ) = ~ J dxe- u ^' D f E dy e u ^l D p Q (y,X). (57) 

— oo 

By multiplying the FPE eq. J5J with e~ At and integrating over time, it is readily verified that po(x, A) appearing 
in eq. (|57|l obeys the following ordinary differential equation 

- 5(x) + Xp = -^(W + D-^)p a . (58) 



Using this equation in the form 
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in eq. I|57|l and integrating a few times by part, the following expression for the first integral in eq. (|57|l is 
obtained 

dye u ^ D p (y, A) = j (o(-x) + D [e u -' D p^ E - e u W D p' Q {x) 

"™'dy U"(y)e u ^/ D p (y) ] j (60) 

where an index "E" means that the respective function is taken at x = xe- For the specific potential eq. I|56(l . 
this yields (using eq. (0) 

dye u M/ D p (y, X) = ^ (&{-x) - e u ^ D p (X) De u ^' D p' Q {x)^ . (61) 

Here, po(\) denotes the Laplace transform of the FPT density for the unperturbed system. Inserting this 
formula into eq. (|52|1 the following quadrature formula is obtained 

h ( XE ) = MR ^e u ^ D p (X) dxe~ u ^l D - J" dxe^W^j . (62) 

One may repeat the whole derivation for arbitrary x (this is needed in the calculation of J2{xe))i yielding 

+ ¥k{f E ' Dpo{X) f d y^ Uiv) ' D ~f dye- u ^/ D Q(-y)\ (63) 

Insertion into eq. I|53[) and a few manipulations of the occurring multiple integrals leads to the following correc- 
tion of the second moment 

j 2 (xe) = -2 + A— 



+ lL£ll poi X)e u ^ D I dx e- u ^/ D I 2 (x) 



where 



X-b dXJ X 

2 / r x E 

X-b 

dx e- u ^' D I 2 {x)-e u ^/ D I{0) dx I{x) ) (64) 

-oo JO ) 



I{x):=e u ^' D dy e- u ^/°. (65) 



B. Formulas for the linear potential case 

I now turn to the specific case of a linear potential, which is particularly simple. Assuming b — and a > 0, I 
have for e — (see, e.g. [37j) 

(T)o = ^f, (66) 
(AT 2 ) = 2^, (67) 



2D 

CV = \j , (68) 

" XeO. 
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and 

/ \\ x E (a- Va 2 +4XD) 
Po(A) = exp — (69) 

I{x) = D/a. (70) 

Inserting the latter expressions into eqs. i|62|) and 164|) yields the following linear corrections of the first and 
second moment, respectively 

h(x E ) = - I ™ - - 1 ] , (71 ) 




2 / X E a D\ x E (a-Va3+4\D) JJ 

32 (x E = s< (-=== + - e 2D _a; B __^. 

' 2 [ Wa 2 + 4AT> a / a J 

(72) 

Hence, the mean and the variance of the first passage time in a linear potential under the influence of a weak 
exponential driving arc given by the fairly simple expressions 

x E (a-^a I +i\D) \ 

■e 2D , (73) 




x E (a-^a 2 +4\D) 
x E )e 2D 

(74) 

It can be easily seen that the corrections to the moments eq. (|71|) and eq. I|72(l are negative if e > (I recall 
that a, D, and A are positive). This makes sense, since a positive force toward threshold will always diminish 
the first passage time and hence also its moments. Furthermore, the correction of the variance is also negative 
for s being positive. 

For very small decay rate (A << 1/(7%), an expansion of eq. l(75|) and eq. (f7^|) in A yields 

(T) -> (r)o-^fA £ , A « 1/(7% (75) 

(AT 2 ) - (AT 2 ) -6^^, A«1/(T) (76) 

This is also obtained if in the formulas of the unperturbed system eq. ()66f) and eq. I|67|) the bias a is replaced 
by a + e\ and the formulas are expanded up to linear order in e. Hence the static limit confirms the result of 
the perturbation calculation. 

Considering the limit of large decay rate (A — > oo), I note that the exponential function in eqs. \T6\ and (|74|l 
can be neglected and I obtain the simple limit 

(T) - A^oo (77) 

(AT 2 ) -> 2D^-^, A^oo (78) 

In the limit A — > oo the driving acts as a 5 spike at t = with amplitude s, which leads to a mod- 
ified initial point x(t = 0) = e. To check this, one can again use the formulas for the unperturbed 
system by replacing the distance between initial point and absorbing boundary (which was x E ) by x E — £■ 
This leads exactly to eq. fTTjl . eq. f7H|l . i.e. in the limit A — > oo the result of the perturbation calculation is exact. 



C. Formulas for the parabolic potential case 



Mean, variance, and Laplace transform of the FPT density for the case of the parabolic potential and e — can 
be written as follows (see, for instance, Refs. pi l37l l4p| 



b ' 



(T) - ^ / dye^erfc(y), (79) 
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(AT 2 ) - ^ / ,///<■">](•(<;)]■-' / ,/: . ■■■(-) (.,- H :) 



(80) 



and 



where 



P_ A/6 (a;_V2j 



a — bxE 



V2Db 



V2Db 



X 2 2 



(81) 



(82) 



In these expressions, erfc(ir) denotes the complementary error function and T> a (z) is the parabolic cylinder 
function |38j . The auxiliary function I(x) reads now 



I(x) 



exp 



(a — xb) 



2 1 



2Db 



erfc [ a , | 
\V2DbJ 



(83) 



Using this function, I obtain the following correction to the mean first passage time out of a parabolic potential 

(84) 



3i{xe) = 



_7T A_ 

2bDX-y 



" x + [e s p (\)erfc(x~) ~ erfc(x+)] 



By means of the derivative of Ji{xe) with respect to A, the correction to the second moment can be brought in 
the following form 



,. i d 1 \ ji(x E ) 2 Xtt 2 

oo 

dx e x2 er{c 2 (x) - e 5 p (X) / dx erfc 2 (x) + 



+crfc(x+) J dx e x erfc(x) 
The resulting formulas for the FPT's mean and variance are given by 



(T) = (T) + 



;A / 7t 



A - 6 V 2bD 



e x + [e p (A)erfc(x_) - erfc(a; + )] 



(85) 



(86) 



(AT 2 ) = (AT 2 ) 



eX 2-k 



X - b V Db 



e 5 erfc(*_) [p' (X) + (T) p (X)} 



+ f dX e^ (xmx - X+ )-e% { X)] 



(87) 



where p' (X) denotes the derivative 01 °f A) (A) with respect to A. 

Again, the limits of small and large-A may be considered. For A — > 0, it is easily seen from eq. (jHJ) that po — > 1 
and po -> -(T) . With this I obtain 



(T) « (T) - ^ ' 



6 V 26L> 



e x -erfc(a;_) - e x +erfc(x + ) , A << 1/(T) 



(88) 



(AT 2 ) m (AT 2 ) 



eXn / 2 



x I I dx e x + x -»*-fr> 2 



erfc z (x) - J dxe x +x +erk z (x) I , A << 1/(T) 



(89) 
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The terms proportional to eX are also obtained by taking the derivative of mean or variance of the unperturbed 
system given in eqs. (|79|) and 1|80|) with respect to parameter a. For very slow driving, the perturbation acts as 
a static change in the bias parameter. Consequently, the perturbation calculation leads to the same result like 
a linearization of mean and variance with respect to a small change in the bias parameter. 
In case of infinite A, the characteristic function po and its derivative p' approach zero yielding the following 
simplified expressions for mean and variance 



(T) - (T) - £ J-^-e<eMx + ) , A — > oo (90) 



(AT 2 ) -> (AT 2 ) - e\j e*+ J dx e* 2 erfc 2 (x), A — > oo (91) 

x + 

Again, what physically happens in this case is a shift of the initial point from x(t = 0) = to x(t = 0) = e 
since the driving e\e~ xt acts as a S spike at the initial time. Consequently, the above results are also obtained 
if in eqs. I|79|l and (|80|l x + (the only term where the initial point enters) is replaced by (a — be)j\/2Db and the 
expressions are expanded to linear order in e. This in turn, is another check that the results achieved cannot 
be completely wrong. 



IV. MEAN AND VARIANCE OF THE FPT: COMPARISON TO SIMULATIONS 



As a verification of the specific results derived in the previous section, I consider the mean, the variance and 
the coefficient of variation (CV) of the FPT in a linear and in a parabolic potential. It will become apparent 
that an exponential driving of these systems can result in a remarkable behavior of the first two cumulants. 
For all data shown, I use a weak positive driving amplitude of e — 0.05, an intermediate noise intensity D = 0.1, 
and the absorbing boundary to be at Xe = 1. Furthermore, two different sets of potential parameters are 
inspected: (i) a = 1, b = for the linear potential and (ii) a = 0.8, b = 1 in case of a parabolic potential. The 
latter choice implies a significantly different FPT statistics since here the escape from the potential minimum 
at x = 0.8 will dominate the passage time (for a large value of a, the potential minimum is beyond the 
absorbing boundary and the FPT statistics will be akin to the linear-potential case) . I compare the analytical 
results derived above to simulation results, that were obtained with a simple Euler procedure in case of a linear 
potential (time step was At = 10~ 4 and 10 6 passage times were simulated), and a modified Euler procedure 
|41j for the parabolic potential (time step was At = 10~ 3 and 10 6 passage times were simulated). In all curves 
data for e — for which I know the exact values of all quantities are shown for the sake of illustration and also 
to demonstrate the validity and accuracy of the numerical simulation procedure. 

Since the agreement between theory and simulations is excellent, I do not have to discuss it at length; I can 
regard it as a satisfying first confirmation of the presented analytical approach. Concerning the agreement 
between theory and simulations at other parameter sets (e.g. larger driving amplitude and smaller or larger 
noise intensity), I restrict myself to the following brief statement: the perturbation result for the correction 
to a moment yields satisfying quantitative agreement with the simulations as long as the correction is small 
compared to the respective unperturbed moment. In general the theory will work best for intermediate up 
to large noise intensities since with a non-weak noise the effect of an additional weak driving will be only 
moderate and the first (linear) correction term will suffice. Note, however, that for systems without a potential 
barrier between initial point and absorbing boundary (like, for instance, the linear potential), the theory works 
at weak noise, too. 

In the remainder of this section I focus on the statistical features of the exponentially driven first passage 
process as they are reflected by mean, variance, and CV as functions of the decay rate A. 



A. Biased random walk with exponential forcing 



It can be expected that a positive exponentially decaying forcing leads to a decrease of the mean FPT. For the 
linear potential (i.e. a biased random walk) this decrease is a monotonous function of A as shown in Fig. ^ 
At small A the correction is proportional to A according to the static approximation eq. i|75|) (shown by the 
dot-dashed line in Fig. . In this range of A the driving is effectively static with amplitude Ae meaning that 
its decay occurs on a time scale that is far beyond the mean FPT. In other words, a part of the driving's power 
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is "wasted" because s(t) still drives the system long after most realizations have been absorbed at x = xe- 
For intermediate values of A, the decay of the driving force takes place much earlier and thus, its effect on 
the mean is less than that of a static driving. Finally, in the large-A limit the mean saturates according to 
eq. (|77|l at the value corresponding to a change of the initial point in the unperturbed system. I note that the 
monotonous behavior of (T) as a function of A differs from what was found for periodically driven linear systems 
in Refs. 0,0. In the latter case, minima @ or maxima Q of the mean FPT vs the driving frequency were 
observed for different boundary and initial conditions. 



1 

<T> 



0.95 
0. 

X 

FIG. 1: The mean of the first passage time vs decay rate of the exponential driving in case of a linear potential. 
Simulations (symbols) and theory (lines) eq. I|73|l with the indicated values of the driving amplitude. The dot-dashed 
line illustrates the slow-driving approximation eq. I75H . 




The variance of the FPT (Fig. |2J) is always below that of the unperturbed case. It shows, remarkably, a 
nonmonotonous behavior as a function of the decay rate. For the parameter set used in Fig. [21 the variance 
attains a minimum at A w 1.6. It is possible to calculate the exact location of this minimum from eq. 174|l and 
express it solely by means of the mean (T)q and the squared coefficient of variation (denoted for brevity by 
R = CVq) of the unperturbed system 



A, ; 



s/l + AR- 2R 2 - 3R 2 /2 +4R-1 
(T) Q R(2~R) 2 



(92) 



0.2^ 



<AT 2 > 
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FIG. 2: The variance of the first passage time vs decay rate of the exponential driving in case of a linear potential. 
Simulations (symbols) and theory eq. (1741 (lines) with the indicated values of the driving amplitude. The dot-dashed 
line illustrates the slow-driving approximation eq. I|7ri[l . 

A minimum in the variance does not occur for an arbitrary parameter set but if and only if 



R < 2 => CV = VWM < v 7 ^, 



(93) 



First passage times under external driving 



14 



i.e. for sufficiently weak noise intensity or large bias a. If the condition eq. (|93|l is met, the value at which the 
minimum is attained is an increasing function of R diverging at R — 2 and saturating for small values of the 
CV (i.e. D -> or a -> oo) at 

^n^lj^r- as CV ->0. (94) 

The occurrence of the minimum seems to be related to the fact that a time-dependent bias reduces the variability 
more strongly than a shift in the initial point (corresponding to the limit A — > oo) does. This gives raise to 
the drop of the variance as A is decreased starting in the large-A limit. The amplitude of the time-dependent 
driving, however, depends on A, too, so its effect on the dynamics gets weaker by further decreasing A and 
the variance starts to increase again. Accordingly, using the exponential driving without the prefactor A (i.e., 
without normalizing the driving's intensity) yields a variance that grows monotonously with A (not shown). 
Therefore, the minimum of the variance is merely based on two competing effects, namely, the greater impact 
of a slow driving (compared to a fast one) on the reduction of the variance and the dependence of the variance 
on the driving's amplitude (i.e. eA). 




X 

FIG. 3: The coefficient of variation of the first passage time vs decay rate of the exponential driving in case of a linear 
potential. Simulations (symbols) and theory using eq. Hl.'-ty ■ eq. I|73jl . and eq. I|74|l (lines) with the indicated values of the 
driving amplitude. The dot-dashed line illustrates the slow-driving approximation using eq. (I13L eq. 1751 . and eq. 1761 . 

The minimum in the variance vs A could be interpreted as an "optimal" decrease in variability due to an 
exponential driving. Things look different, though, from the view point of relative variability as it is quantified 
by the CV (cf. Fig. |3J). First of all, depending on the value of A, the CV can be both larger or smaller than 
in the unperturbed case. For the value of A at which the variance attains a minimum the CV is larger than in 
the unperturbed case and is thus far from being "optimal" . The fact that the effect of the exponential driving 
on the CV can be both positive or negative can be understood by looking at the CV in the limit cases of small 
and large decay rate where it corresponds to the CV of the unperturbed system with rescaled parameters. The 
latter depends on the inverse of the product axE- A static increase of the bias (replacing a by a + eX which 
is the effect of a slow driving) will thus lead to a decrease in CV compared to the unperturbed case (cf. the 
static approximation in Fig. In contrast, diminishing the difference between initial point and absorbing 
boundary (replacing xe by xe — £ which is the effect of the forcing for A — > oo) leads to a higher CV than in 
the unperturbed case. Interpolating between the two limit cases will inevitably lead to at least one minimum 
of the CV vs A. For the parameter set used in Fig. [3J this minimum is attained at a decay rate that is smaller 
than the inverse mean FPT of the unperturbed system. 

In conclusion, already in the simple linear potential case, the effect of an exponential driving can be fairly 
involved. For the behavior of mean, variance, and CV as functions of the decay rate, it was essential that I have 
used a constant-intensity scaling of the driving function. 

B. Escape out of a parabolic potential with exponential forcing 

With b > there is a true state dependence on the right hand side of the dynamics eq. The state variable 
is attracted toward the potential minimum at x — a/b. If this rest state is far beyond the absorbing boundary 
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xe (i-e. xe <§C a/b), the FPT statistics will be similar to that of the biased random walk. In the following 
I choose, however, b = 1, a = 0.8, and x E — 1 such that a/b < xe- With this choice the FPT problem 
is significantly different from the linear potential previously discussed, since in order to reach the absorbing 
boundary at xe — 1 the state variable x(t) has to be driven by a sufficient noise to accomplish the escape out 
of the potential minimum. 

In case of a parabolic potential, already the mean FPT depends nonmonotonously on the driving's decay rate 

2.7, 1 




X 

FIG. 4: The mean of the first passage time vs decay rate of the exponential driving in case of a parabolic potential. 
Simulations (symbols) and theory eq. (1861 (lines) with the indicated values of the driving amplitude. The dot-dashed 
line illustrates the slow-driving approximation eq. (18811 : the dotted line is the extended large- A approximation eq. 1981 . 



A (Fig. |3J; it attains a minimum for A w 1.5 which stands in marked contrast to the linear case. The reason 
for the occurrence of this minimum is the state dependence of the dynamics eq. QJ as I will show now. First of 
all, starting at A = 0, the mean FPT decreases linearly with A in accord with the static approximation shown 
by the dot-dashed line in Fig. 0] This is completely equivalent to the linear potential case. Second, the mean 
FPT also drops if the decay rate is decreased starting from the large-A limit. In other words, the limiting value 
is approached from below. The behavior of the mean in these two limits implies the occurrence of at least one 
minimum. 

In order to understand why the large-A limit is approached from below, I consider a large but finite value of A 
and a time t for which 

j « t « (T) . (95) 

For such times it can be taken for granted that the driving has practically decayed to zero but on the other 
hand it is highly unlikely that a realization has already reached xe ■ In this case the effect of the driving can be 
inferred from the free solution of eq. Q for a parabolic potential with initial value x(t — 0) — which can be 
written as follows 

t 

x(t) =e(l+ j^-^j e- bt - £— ^e~ At + J dte b ^ [a + VziJf (i)] (96) 

o 

If the condition eq. (|95|l is met, the second term can be neglected and the solution reads 

t 

x(t) =s(l+ e- bt + J dte b{i - t] [a + y/2D£(t)] (97) 

o 

This approximate solution is, however, equivalent to the unperturbed dynamics with an initial point at x(t — 
0) = e[l + b/(X — b)]. The equivalence holds true for a time t obeying eq. (|93)l and any time larger than this 
time. In other words, for t > 1/A the realization of the original process and that of the unperturbed process 
with the modified initial condition differ only by a small exponential contribution. Consequently, also the FPT 
statistics of both processes will be the same provided that a successful escape toward xe is highly unlikely for 



First passage times under external driving 



16 



short times at which eq. I|97|) does not hold true. 

For A — * oo, the initial condition approaches the value x(t = 0) — > e (see the discussion around eq. (|§UJ)). For 
a large but finite value of A, the shift in the initial point will be larger than in the latter limit and thus the 
mean will be more strongly decreased than in the strict limit A — > oo. To obtain an explicit formula showing 
this drop, the mean of the unperturbed system with modified initial point is linearized with respect to e (this 
is not strictly necessary but consistent with the linear approach used throughout this paper) yielding 

(T)«(T) -e(l + x ^)y^e*+erfc( a:+ ), A » 1/(T) (98) 

Of course, this is also obtained by replacing the amplitude s in eq. (|90[l by the modified amplitude e[l+b/ (A — b)]. 
The approximation is shown in Fig. 01 by a dotted line; it displays the drop of the mean with decreasing A at 
large decay rate and agrees well with the full solution in this range. 

I note that a nonmonotonous behavior of the mean was also found for a system with parabolic potential and 
a periodic forcing in ref. |15| . If the system is driven by s(t) = £COs(uit) the mean passes through maxima and 
- less pronounced - minima when plotted as a function of the driving frequency (cf. in particular Figs. 11,12, 
and 13 in [T^V There are two important differences between the exponential and periodic driving functions: 
(1) the periodic driving attains both positive and negative values; (2) the amplitude of the periodic forcing as 
considered in Ref. |15| is fixed and does not depend on the time scale of the driving. The maxima and minima 
found for periodic driving are true resonance phenomena. In contrast, the minimum in the mean FPT for 
exponential driving appears as a compromise between the dependence of the driving's amplitude on the decay 
rate (drop of the mean with increasing A at small A) and the stronger effect of a truly time dependent driving 
on the state-dependent dynamics (drop of the mean with decreasing A for A 3> 1/(T)q). I would like to point 
out that the latter argument does not apply in case of a linear potential (i.e., a state-independent force) because 
for b — the A dependent modification of the initial point in eq. I|97[) vanishes. Hence, in this case we cannot 
infer the existence of a minimum of (T) vs A and, in fact, it also does not occur as was seen in the previous 
subsection. 
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FIG. 5: The variance of the first passage time vs decay rate of the exponential driving in case of a parabolic potential. 
Simulations (symbols) and theory eq. (1871 (lines) with the indicated values of the driving amplitude. The dot-dashed 
line illustrates the slow-driving approximation eq. (18911 : the dotted line is the extended large- A approximation eq. 1991 . 




Turning to the variance depicted in Fig. I note that this function also passes through a minimum vs A like 
in the linear case. This minimum occurs at a smaller decay rate (A « 0.38) than that of the mean FPT and 
remarkably close to the inverse mean first passage time of the unperturbed system (1/(T) « 0.37). Plotting 
the analytical solution eq. (|87() for different parameters revealed that this time-scale matching condition holds 
true as long as the system is in the "subthreshold" regime, i.e. for a/b < xe and weak up to moderate noise 
intensity. For larger noise intensity and/or "suprathreshold" system parameters (a/b > xe) the minimum is 
attained at values larger than l/(T)o. For the specific limit of weak noise and a — > oo, one can expect the 
location at the value found for the linear potential, namely \ m in — ► 3/(2(T)o) which is indeed larger than the 
inverse of the mean FPT in the unperturbed case. 

The minimum of the variance can be understood by the same line of reasoning as in case of the mean, i.e. by 
considering the behavior at small and large-A which are determined by the static approximation and by the 
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effective soiution eq. (|97|) . respectiveiy. the latter leads - in complete analogy to the derivation of eq. I|98[l - to 
the extended large- A approximation of the variance 

. oo 

(AT 2 ) « <AT 2 ) -e(l+ j^j yj j^*e*+ J dx e^erfc 2 ^), A » -L (99) 

x + 

This is shown in Fig.|53by the dotted line. I note, however, that the actual drop in variance at large A extends 
over a much larger range where eq. I|99|l does not hold true anymore; the decrease of the variance in this range 
is also much stronger than expected from eq. Ij99(l . The effect of a temporally extended driving is thus much 
stronger than a change of the initial point similarly to the case of a linear potential. Furthermore, because the 
amplitude of the driving depends on A, the variance will drop for A — » to the value of the unperturbed system. 
The occurrence of the minimum is therefore mainly based on the different sensitivity of the FPT statistics with 
respect to changes in the bias term a and the initial point of the passage and the A dependence of the driving 
amplitude due to the constant-intensity scaling of the forcing. 



0.69 
0.685 

cv 

0.68 

0.675 

0.67 
0. 

X 

FIG. 6: The coefficient of variation of the first passage time vs decay rate of the exponential driving in case of a parabolic 
potential. Simulations (symbols) and theory using eq. (1131 . eq. 18611 . and eq. 18711 (lines) with the indicated values of the 
driving amplitude. The dot-dashed line illustrates the slow-driving approximation using eq. 113L eq. 1881 . and eq. 1891 1. 




Since the minima in mean and variance vs A are attained at distinct values of the decay rate, I can expect a 
complicated behavior for the relative variability of the FPT. Indeed, the CV as a function of decay rate shown 
in Fig. first goes through a minimum (X m in ~ 0.66), reaches a maximum at a finite decay rate (X max ~ 2.8), 
and saturates at a CV that is higher than in the unperturbed case. Compared to the CV of the linear potential 
system (cf. Fig. I3J), the decrease at small A is weaker; furthermore, there is no maximum for the linear system 
but only a saturation at large decay rate. 

For the chosen parameters the FPT from x = 0tox — xe can be split into two independent FPTs as T — T\ +T2 
with T\ being the FPT from x — into the minimum x = a/b and T2 being the time for the passage from the 
minimum to the absorbing boundary (see, for instance, ref. 142(1. It is straightforward to show that 

cy2 = cv ^ + cv^ (100) 

where CV\ and CV2 are the CV of the respective passage processes. Now the relaxation into the minimum 
is evidently more regular than the noise-assisted escape out of the potential minimum, i.e. CV2 > CV\. The 
behavior of the CV can be understood by considering how in the limits of small and large decay rate the relative 
contributions of Ti and Ti are changed. 

At low decay rate, the driving is effectively static, hence the system is equivalent to the unperturbed dynamics 
with enlarged bias a + eX. This system in turn is equivalent to the unperturbed dynamics with the original bias 
a but initial point at x = — eX and absorbing boundary at x — xe — sX. With these parameters, the FPT T\ 
from initial point to potential minimum increases and the FPT T2 from minimum to absorbing boundary drops 
compared to the unperturbed case. It is reasonable that the CV of the single passage processes change only 
little; what mainly changes is their relative contribution to the total CV by means of the squared ratios T\/T 



First passage times under external driving 



18 



and T2/T. Hence, according to eq. (|100fl the CV will drop since the lower CV of T\ makes a larger relative 
contribution. From this line of argument, it is also reasonable that the CV drops with increasing A as long as 
the static approximation holds true. 

In the limit of A ^> (T), one obtains also the unperturbed dynamics (as explained above by means of eq. i|9 1\\ 
but for a passage process from x = e[1 + b/(X — b)] to x = xe- Obviously, here the escape time T2 is the same as 
in the unperturbed case; the first FPT Ti, however, has been shortened. According to eq. \ 100(1 . one can thus 
expect a higher CV than in the unperturbed case since CV2 makes a larger relative contribution to the CV of 
the total FPT. Moreover, the A — > 00 limit of the CV is approached from above because the shift in initial point 
drops with increasing A. Interpolating in the simplest way between the behavior at small and large A predicts a 
minimum at moderately low decay rate and a maximum at larger rate as has been found in Fig. I note that 
for the decay rate at which the variance attains a minimum, the CV is higher than in the unperturbed case, 
similarly to the linear potential case (b = 0) discussed in the previous subsection. 

Like in case of the biased random walk, the behavior of mean and variance as functions of the decay rate depends 
crucially on the constant-intensity scaling of the driving I have used. Additionally, the state-dependence of 
the force leads to a nonmonotonous behavior of the mean as a function of the decay rate. The minima in 
mean and variance are not true resonances as in case of a periodic driving but are mainly related to the 
distinct sensitivity of the FPT statistics with respect to changes in the initial point or in the bias parameter, 
respectively. Nevertheless, in physical situations where a constant-intensity scaling of the driving is appropriate, 
the nonmonotonous behavior of the first two cumulants and of the CV may be of some importance. 

V. SUMMARY AND OUTLOOK 

In this paper I have studied the moments of the first passage time in presence of a weak time-dependent driving. 
A formula for the corrections to the moments for arbitrary driving and potential shape were derived that 
contains, however, the time-dependent probability density of the unperturbed system or its Laplace transform. 
The latter functions are known only in a few rare cases. Explicite correction formulas for mean and variance of 
the first passage time could be achieved for the case of an exponentially decaying driving function and an either 
linear or parabolic potential. These analytical results were found to be in excellent agreement with results from 
computer simulations of the passage processes. I demonstrated furthermore, that for the chosen exponential 
driving, the variance of the passage time in the linear case as well as both the mean and the passage time in 
case of a parabolic potential pass through minima as functions of the decay rate of the driving. The behavior 
of the relative standard deviation (i.e. the CV) proved to be even more complicated. All of these findings 
resemble the effects of coherent stochastic resonance and resonant activation, however, they rely on a different 
mechanism involving the constant-intensity scaling of the driving. 

The explicit results for the case of a linear or a parabolic potential derived in this paper will be applied in the 
near future to the neurobiological problems mentioned in the introduction. Furthermore, the results can be 
useful for other problems, too. The application of the formulas to the case of a periodic driving with a cosine 
function studied in Refs. 0, 0, is straightforward. Likewise, the case of an exponentially damped cosine 
function involving two time scales (driving period and decay rate) can be readily derived from my formulas and 
might be worth to look at. 

The general approach presented in this paper may be easily extended to the cases of two absorbing boundaries 
or one absorbing and one reflecting boundary. I am also convinced that the problem of a state-dependent driving 
(i.e. dealing with a force s(x, t) instead of s(t)) can be successfully treated with the approach. Finally, the case 
of a stochastic driving function might be tackled by a proper average of the correction formulas over the driving 
process and its initial condition. This last problem, though, seems to be much more challenging than the other 
extensions of the theory. 
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APPENDIX A: EQUIVALENCE OF THE DIFFERENT QUADRATURES EXPRESSIONS 

Here I show that eq. (|27f) together with eq. (|24ll and eq. I|16f) yields the same FPT moments as the standard 
formula eq. 10. The two differing expressions for the n-th moment can be written as follows (indices "S" and 
"A" stand for "standard" and "alternative") 

b yi b V2 

( T )s = d yi e d dxie dy 2 e ° dx 2 e ° ■■■ 

a — oo xi — oo 

b y n 



dy n e~^~ / dx n e ~ (Al) 



x n -l 



b b x\ b 

(T n ) A = — I dy x e ° dx x e dy 2 e ° dx 2 e b 



yi -oo 

b 



D n 



dy n e d Jdx n e ^Q(x n -a) (A2) 

y n 

Note that for the second formula, I have used an arbitrary initial point x(t = 0) = a which only changes the 
argument of the Heaviside function in eq. (|24|) . 
Now I introduce the operators 

b y 

t Xj {x) = fdye u ^/ D f dxe~ u ^' D (A3) 



-OO 

b 



M x Ax) = / dy e - u ^/ D / dxe u ^' D (A4) 



where the argument indicates the variable with respect to which the respective function is integrated, while the 
index denotes as a parameter one boundary of integration. It is not hard to show that for u,w < b 

M u (v)Q(v-w) = t w (v)S(u-v). (A5) 

and 

M u {v)t w {x) = t w {x)M u {v), (A6) 
i.e., the operators commute if their arguments and indices differ. 

By means of the operators, the two expressions for the n-th moment can be written as follows 

= Tsr /C«(a:i)/C a!l (a:2)---/C 1BB _ 1 (i„)e(6-a: n ) (AT) 
< T ">^ = -^M b {x x )M xl {x 2 )---M Xn _,(x n )Q{x n -a) (A8) 

L) 

Note that the function 9(6 — x n ) in eq. I|A7|I is always one since Xj < b for all j = 1, . . ,n. Using the relations 
eq. (|A5|) and eq. (|A6|> it follows from eq. (|A8|) that 
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